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Abstract. Lately, many phenomena in both apphed and abstract mathemat- 
ics and related disciplines have been expressed in terms of high order and frac- 
tional PDEs. Recently, AUouba introduced the Brownian-time Brownian sheet 
(BTBS) and connected it to a new system of fourth order interacting PDEs. 
The interaction in this multiparameter BTBS-PDEs connection is novel, leads 
to an intimately-connected linear system variant of the celebrated Kuramoto- 
Sivashinsky PDE, and is not shared with its one-time-parameter counterpart. 
It also means that these PDEs systems are to be solved for a family of func- 
tions, a feature exhibited in well known fluids dynamics models. On the other 
hand, the memory-preserving interaction between the PDE solution and the 
initial data is common to both the single and the multi parameter Brownian- 
time PDEs. Here, we introduce a new — even in the one parameter case — 
proof that combines stochastic analysis with analysis and fractional calculus 
to simultaneously link BTBS to a new system of temporally half-derivative 
interacting PDEs as well as to the fourth order system proved earlier and dif- 
ferently by AUouba. We then introduce a general class of random fields we 
call inversc-stablc-Levy-time Brownian sheets (ISLTBSs), and we link them to 
/3-fractional-time-derivative systems of interacting PDEs for < /3 < 1. When 
/3 = l/u, V S {2,3,...}, our proof also connects an ISLTBS to a system of 
memory-preserving f-Laplacian interacting PDEs. Memory is expressed via 
a sum of temporally-scaled fc-Laplacians of the initial data, fc = l,...,i/ — 1. 
Using a Fourier-Laplace-transform-fractional-calculus approach, wc give a con- 
ditional equivalence result that gives a necessary and sufficient condition for 
the equivalence between the fractional and the high order systems. In the one 
parameter case this condition automatically holds. 



Date: January 20, 2013. 

1991 Mathematics Subject Classification. Primary 35C15, 35G31, 35G46, 60H30, 60G60, 
60J45, 60J35; Secondary 60J60, 60J65. 

Key words and phrases. Brownian-time Brownian sheet, linear systems of fourth order inter- 
acting PDEs, linear systems of high-order interacting PDEs, linear systems of fractional interacting 
PDEs, Brownian-time processes, initially perturbed fourth order PDEs, memory-preserving PDEs, 
Brownian-time Feynman-Kac formula, iterated Brownian sheet, random fields. 

1 



2 



HASSAN ALLOUBA AND ERKAN NANE 



Contents 

0. Introduction and statement of results [3 
0.1. Recalling the Brownian-time Brownian sheet and its fourth order 

interacting PDEs system 3 

0.2. The new BTBS time-fractional interacting PDEs system |5 
0.3. Conditional equivalence of the BTBS fourth order and time-fractional 

interacting PDEs systems @ 
0.4. Inverse-stable-Levy-time Brownian sheets and their interacting 

fractional and high-order PDEs systems 
0.5. Conditional equivalence of the ISLTBS high order and time-fractional 

interacting PDEs systems IC 

0. 6. Examples of the initial function / 11 

1. Establishing the BTBS fourth-order and fractional PDEs systems links 12 

1.1. A stochastic analytic unifying proof of Theorem 10. II and Theorem 10.01 12 

1.2. Proof of the conditional equivalence result: the BTBS case 16 

2. Establishing the ISLTBS high-order and fractional PDEs systems links 18 

2.1. A stochastic analytic proof of Theorem 10. 31 18 

2.2. Proof of the conditional equivalence result: the ISLTBS case 23 

3. On relaxing the boundedncss condition on / and its derivatives 25 
Acknowledgement [26 
Appendix A. Another proof of the fractional PDEs connection in 

Theorem [03] [26 

Appendix B. Frequent acronyms and notations key 27 

References 27 



0. Introduction and statement of results 

Many phenomena in mathematical physics, fluids dynamics and turbulence mod- 
els, mathematical finance, and the modern theory of stochastic processes were re- 
cently described through fractional and higher order evolution equations (e.g., see 
[1]-[TT], [13], [l5j-[18], and [20]-[34]). Another feature important in fluids dynamics 
and other applied fields is the interaction between different quantities of interest 
(velocity and pressure in the Navier-Stokes example). In his article [1], AUouba 
introduced n-parameter d-dimensional Brownian-time Brownian sheet (BTBS) — a 
Brownian sheet where each "time" parameter is replaced with the modulus of an 
independent Brownian motion — and connected it to a new system of fourth or- 
der interacting PDEs. The interaction in these PDEs added novel and intricate 
new features not present in their one parameter {n ~ 1) c?-dimensional {d > 1) 
Brownian-timc-Brownian- motion PDE counterpart — first given in [6l|5] (DeBlassie 
also gave a different proof in [18]). This interaction means, among other things, 
that the PDEs are to be solved for a family of ti + 1 functions and not a single func- 
tion as in the one-time-parameter case. As in the n ~ 1 case, the BTBS-PDEs are 
memory preserving. I.e., the solution family for the BTBS fourth order PDEs also 
interacts with the initial data via a time-scaled Laplacian of the initial function. 

In this article, wc first introduce a new — even in the one parameter case — proof 
that judiciously combines Ito's rule with properties of the Brownian sheet as well 
as analysis and fractional calculus to simultaneously link BTBS to a new system of 
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temporally half-derivative interacting PDEs as well as to the fourth order system 
proved in [Tj by Allouba via a different approach. Like the fourth order PDEs 
connection, the interacting feature of this half-derivative PDEs system belongs 
solely to the multiparameter setting and docs not exist in its one-parameter coun- 
terpart. Here, we note that the half-derivative connection to the one-parameter 
version of BTBS — the Brownian-time Brownian motion — was first noted and es- 
tablished by Allouba and Zheng in [6] through their half-derivative generator, an 
implicit equivalence of the Brownian-time Brownian motion fourth order PDE to 
a half-derivative-in-time fractional PDE. The fractional PDE connection and its 
equivalence to the fourth order PDEs was then given explicitly and generalized by 
Nane, Meerchaert, Baeumer, Vellaisamy, Orsingher, and Beghin [ST j [29 l [T5 l l32 | [T6 ] . 
In the second part of the paper, we introduce the class of inverse-stable-Levy-time 
Brownian sheets (ISLTBS's) and, using an adaptation of our proof for the BTBS, we 
link any given ISLTBS to a /3-fractional-time-dcrivative system of interacting PDEs 
for < /3 < 1. When, (3 = l/v, z/ e {2,3, . . .}, our proof also connect our ISLTBS 
to a system of high order interacting PDEs involving the spatial i^-Laplacian tcrn|3 
I^^. Like the BTBS case, the interacting nature of these PDEs is unique to the 
multiparameter setup and it vanishes in the one-time parameter framework. The 
memory-preserving feature for ISLTBS is apparent in its high (2;^) order PDEs via 
a sum of temporally-scaled spatial fc-Laplacians of the initial data, k = 1, . . . ,v — 1. 
In all cases, using a Fourier-Laplace-transform-fractional-calculus approach, we give 
a conditional equivalence result that gives a necessary and sufficient condition for 
the equivalence between the fractional and the high order systems. In the one pa- 
rameter case this condition is automatically satisfied. In a separate article |10| . we 
treat the case of stablc-Lcvy-time fractional Brownian sheet connection to PDEs 
systems. 

We now introduce the setup and state our main results. It is both instructive 
and motivating to first introduce and treat the Brownian-time Brownian sheet case. 

0.1. Recalling the Brownian-time Brownian sheet and its fourth order 
interacting PDEs system. Let _B'^' , . . . , i?^") be n independent copies of a 
standard one-dimensional Brownian motion starting at and independent of an 
n-parameter (n > 1) M''-valued Brownian sheet 



1^0 ^ {lyO(t) = (W^O(t), . . . , I¥rf"(t)) ;t ^ (ti, . . . ,t„) e M!^} , 

"starting" at e R*^ under P— P{W^''(t) = 0} = 1 for t e dW% (see Notation [Ol])— 
all defined on a probability space (f2,^,P). Of course, the Brownian sheet co- 
ordinates {VF°(t); t > O} , . . . , {iV^(t); t > O} arc assumed independent. For any 
x={xi,...,Xd) £ M'', let 

(0.1) W^^ ■.= W'' + x=^{W^{t) = (M^fnt),---,W^?(t));t = (ii,...,i„) €M!^}. 
The Brownian sheet transition density is given by 



(0.2) 




Frequently used notations and acronyms can bo found in Appendix iBl 
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We define the 71-parameter R'^-valucd Brownian-time Brownian sheet (BTBS) 
based on and B'^^\ . . . and starting at a; e M^, by 



(i?W(ti) 



t e 



(0-3) W^,„ (t):= 

Clearly, P W^(i) g(„)(t) = x = 1 for t e (9M". Our first main result gives a 
nonlinear fourth order interacting PDE and a corresponding linear system of fourth 
order interacting PDEs that are solved by running the BTBS in (jO.Sp . Before stating 
our main result, it is helpful to adopt some simplifying notational conventions. 

Notation 0.1. We always use the notation N„ :— {!,... ,11}; and we denote by 

o 

M" and the interior and boundary of R" , respectively. We will alternate 

freely between ut, it(t), and u(ti,...,t„) for typsetting convenience and ease of 
exposition. Moreover, ty= (ti)igN^^\{j} S R"""'^ will denote the (n — l)-tuplc point 
in M""^ that is obtained from t = (ii, . . . , t„) S R" by removing the j-th variable, 
tj. The notations Ut, it(t), and uty{tj) will all mean the n-parameter function 
u : — >• E** evaluated at t = (ti, . . . ,i„); and we use uty{tj) whenever we need 
to focus on the changes in u as a function of the j-th variable tj while holding 
the rest of the parameters in ty fixed (e.g., when we apply Ito's rule in tj). The 
same comment applies for the notations u{t,x) and uty{tj,x). In denoting partial 
derivatives of any positive integer order k in the variable i^, we use d^. and the 
mixed derivative in ti, . . . , i„ by . 

We recall first AUouba's fourth order interacting BTBS-PDEs connection. Here, 
we only state the linear system connection (see [1] for the correspondig nonlinear 
fourth order interacting PDE). 



Theorem 0.0 (AUouba [I]: BTBS and its fourth order interacting PDEs con- 
nections). Let {W^(i) g(„)(t);t g R"} he an n-parameter -valued BTBS 

based on a Brownian sheet and Brownian motions {^B^^^^ and starting 
on dM.\ at X on (17, ^,P). Let f -.W^ he hounded and measurahle 

such that all second order partial derivatives d^^ ^^f are bounded and Holder 
continuous, with exponent < a < 1, for 1 < k,l < d. If 



w(t,x) ==Ep [/ (w^a,,...,B(„,(t) 



(0.4) 



'^(J)(t,a;) = 



v«eN„\{j} 



(t) 



for (t, x) e R![ X R'^, then the family {u, '^'^^^ } 
of interacting fourth order linear PDEs: 



is a solution to the system 



(0.5) 



(a) dt^u 



»ew„\{j} 



A^f + -Al9/^^^; te R" ,xe 



(b) u{t,x) = f{x); 

(c) '^''^^t,x) = 0; 

(d) ^(j)(t,x) = [n. 



1)71—1 



t e ( 

tj = 0,.T e 



, X e 
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forj e N„. 

The intriguing interaction (or coupling) in Theorem 10.01 between u and 
j ~ 1, . . . ,n, shows that the PDEs in (|0.5p are nontrivial intricate generahzations 
of the Brownian-time Brownian motion PDE connection in [6j[5] (the case n = 1), 
given by 

[u{0,x) = fix); cceM'^, 

in which u ~ This couphng phenomenon in the rt > 1 is caused by the 

interaction between the Brownian-times jiJ^^'' | , . . . , li?^"' | and the variance of the 
outer Brownian sheet through c^t^ Var (ly^(f)) = nieN„\{j}*'- '^^^^ its 
one-parameter Brownian-time process (BTP) counterpart, the BTBS is memory 
preserving as is indicated by the inclusion of the Laplacian of the initial function 
A^f{x) in the BTBS PDEs in dHS]). This means that the PDE's dynamics in (|(l5)) 
depend on f{x) for tj > 0, and the role of / carries over beyond its more customary 
starting time — tj = — in more traditional PDEs. This preservation of the effect of 
the initial data f{x) is a manifestation of the non-Markovian nature of the BTBS 
in the time parameter tj. The effect of the initial data, through A.xf, tapers off as 



tj — 7> oo at the rate of y (JlieN "tjTr") . Also, just as BTPs are not 

classical (not semimartingales, not Markovian, and not Gaussian), BTBS is not a 
classical random field. 

0.2. The new BTBS time-fractional interacting PDEs system. The Caputo 
fractional derivative of order < /? < 1, defined by 

(0.7) d^u{t,x) = dMs,-)j^, 

whenever u is in the time parameter. One of the main advantages of the Caputo 
fractional derivative is that initial conditions are expressed in terms of initial values 
of integer order derivatives. Thus, the Caputo fractional derivative is well suited for 
the Laplace transform techniques and to properly handle initial values [T7]. It has 
been widely used to solve ordinary differential equations that involve a fractional 
time derivative [22l |34] . 

Our first main result for the BTBS provides its new link to a temporal half- 
derivative PDEs system. 



Theorem 0.1 (BTBS and its interacting fractional PDEs system). Assume 
that {W^(i, ^(„)(t);t S M"} is an n-parameter W^-valued BTBS based on 

a Brownian sheet and Brownian motions {_B('-'}"_^ and starting on 9M" 
at X G M.'^ on (fi,^,P). Assume the same conditions of Theorem \0.0\ on 
f :R'^ ^R. Ifu is as in TheoremWll and 



(0.8) r(^')(t,.T) =Ef 



k*6N„\{j} 
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for (t, x) e R" X M'^, then the family {m, '^(•J)} is a solution to the system 



(0.9) 



(a) a.^.ii- ^A.rW; 



t e R" , X e 



t G 



a; e 



1 

(fe) M(t,a:) =/(.t) 

(c) ^^(t,^) =0; ^ tyG...+ 

jd) r«(t,z) = [u^eN^\{J} - 0,.T e 

for j e N„. /n particular, if n = 1, u = ^(j) and (|0.8p reduces to 

i (a) d^u = ^A.^u; t £ (0, oo), x e R^, 
(0.10) J V / t ^ ^ ' ^ ■ ^■ 

I (6) u(0,x) ==/(x); xeR^ 



x G R"*, 



Our proof of this result is stochastic analytic in flavor, and it is illuminating since 
it unifies the fourth order and fractional PDEs in one new equation (see Subsection 
12. II below). In doing so, it 

(i) simultaneously establishes Theorem 10.01 and Theorem 10.11 for the BTBS, 
and in the process it gives another proof to Theorem 10.01 different from 
that given in [T] (this also provides a novel unifying way to concurrently 
prove the fourth order and fractional PDE connections to the one parameter 
(n — 1) Brownian-time Brownian motion case); and it 

(ii) leads to a revealing proof (new even in the one parameter n ~ 1 case) in 
the more general case of inverse Levy-time Brownian sheet case that links 
the ISLTBS to two new interacting PDEs systems: one fractional in time 
and the other of order 2i/, ly e {2, 3,4,.. .}. 

It is now important to highlight the significant difference between the fourth- 
ordcr-fractional PDEs connection in the one-parameter Brownian-time Brownian 
motion case and that in the multiparameter Brownian-time Brownian sheet case 
here. Unlike (|0.6p and (|0.10p . which are both single PDEs in one function u, 
the n > 1 cases in ()0.5p and (|0.9p are systems of interacting PDEs involving two 
different families of functions |u, ^-' j^gj^ and {u, ■ This means that an 

equivalence between the fourth order system (jO.Sp and the half-derivative system 
([O)) may only be established by relating and Y^^'^ for j G N„. This is the 
content of our second main result given in the following subsection. 

The Caputo half-derivative feature of absorbing the initial condition u{0,x) in 
PDEs as compared to the first derivative (compare (j0.5|) with (|0.9|) ) can be seen in 
several different ways. The first way is through our stochastic analytic proof which 
simultaneously proves Theorem 10.01 and Theorem 10.11 The second way is through 



its Laplace transform s^u{s,x) 
as well. 



s^ ^u{0,x) which incorporates the initial value 



0.3. Conditional equivalence of the BTBS fourth order and time-fractional 
interacting PDEs systems. Our second main theorem for the BTBS furnishes 
the condition needed for equivalence to hold between the fourth order and the 
fractional systems in Theorem 10.01 and Theorem 10.11 
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Theorem 0.2 (Conditional equivalence of the fractional and the fourth-order 
interacting PDEs systems). Assume that f satisfies the conditions of Theo- 
rem 1 0. 0\ and that the pair , ^(^^} satisfies 



' (a) %/8A,(ajy'r(^')) = A2^(^) ; t £ M'^ , a; £ 



(0.11) 



(c) '2^(j')(t,a;) = 
(e) r(^')(t,.T) = 



n?eN„\{j} 

^(d) r^3\t,x) = '?/(^')(t,.T) =0; tye d^X^^x e M'', 

/or some j* G N„ ('a/Z t/ie derivatives exist and are continuous). Then (|0.9p 
is satisfied for j* by {u^Y^^ '} ijff (|0.5p is satisfied for j* by {u, ^'-^ 

Conversely, if (jO.Sp and (|0.9p are satisfied by 

/or J* £ N„, t/ien (jO.lip /loMs for j* . In particular, the fourth order 
system { (|0.5p }^^f,j and the fractional system { (j0.9p }j^f.j are equivalent iff 
the collection of functions {'^■'V^^-'' j^gp^ satisfies (jO.lip /or j G N„. // 
{'^■' , l^^^-* are t/ie BTBS functionals defined as in Theorem 1 0. 6^1 and T/ie- 

orem W.li then {'^\y^^'>} .^^ satisfy (|0.1ip /or every j G N„. Ifn = l, then 

one-parameter version of the condition (jO.lip is satisfied, 
and the equivalence between (j0.6p and (jO.lOp holds. 



0.4. Inverse-stable-Levy-time Brownian sheets and their interacting frac- 
tional and high-order PDEs systems. Inverse stable subordinators — which we 
also call inverse Levy motions — arise in |281 127j as scaling limits of continuous time 
random walks. Let S{n) ~ Yi + ■ ■ ■ + Yn a sum of independent and identically 
distributed random variables with EYn = and EY^ < oo. The scaling limit 
c~^^'^S{[ct]) B{t) as c — > oo is a Brownian motion B at time t, which is normal 
with mean zero and variance proportional to t. Consider y„ to be the random jumps 
of a particle. If we impose a random waiting time r„ before the nth jump Yn, then 
the position of the particle at time T„ = Ji + • • • + J„ is given by S{n). The number 
of jumps by time t > is N{t) ~ max{n : T„ < t}, so the position of the particle 
at time i > is S{N{t)), a subordinated process. If P(J„ > t) = t~^l{t) for some 
< /3 < 1, where l{t) is slowly varying, then the scaling limit c^^/'^Tjctj => L(t) is 
a strictly increasing stable Levy motion L at time t and with index /3, sometimes 
called a stable subordinator. The jump times r„ and the number of jumps N(t) are 
inverses {N{t) > x} = {T([x]) < t} where [x] is the smallest integer greater than 
or equal to x. It follows that the scaling limits are also inverses c~^N[ct) ^ A(t) 
where A(t) = inf{x : L{x) > t}, so that {A{t) < x} = {L{x) > t}. We call the pro- 
cess A a /3- inverse Levy motion. Since N{ct) k, K{t), the particle location may, 
for large c, be approximated by c-'^/'^S{N{[ct])) w {c'^)~^/'^ S(ct^ K{t)) sa B{A{t)), a 
Brownian motion subordinated to the inverse or hitting time (or first passage time) 
process of the stable subordinator L. The random variable L{t) has a smooth den- 
sity. For properly scaled waiting times, the density of L{t) has Laplace transform 
e"**''' for any t > 0, and the random variables L{t) and t^^^L{l) are identically 
distributed. Writing gp{u) for the density of L{1), it follows that L{t) has density 
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t-'^/f^gpit-^/'^u) for any t > 0. Using the inverse relation P(A(i) < a;) = F{L{x) > t) 
and taking derivatives, it follows that A{t) has density 

(0.12) K^/^^ = tr'x-'-^/'g,{tx-'/^), 

whose i i~> s Laplace transform s^~'^e~^^'^ can also be derived from the equation 

by taking Laplace transforms on both sides. Some fundamental properties of K^^f^ 
are summarized below 

Lemma 0.1 (Hahn ct al. [13). Let K^:^^ be the density of K{t). Then 

(a) \\mt~^j^Q K^.^^^ = 5q{x) in the sense of the topology of the space of tempered 
distributions ^'(K); 

(b) \\ui^^+qK^£^ = r^-^, t > 0; 

(c) Um.^oo = 0, i > 0; 

(d) t I— > s Laplace transform of K^'q^ is s^~^e~^'^^ ; 

(e) For each t > 0, K^^^^ satisfies 

(0.13) <oi - -a.i^X - Y^^^^oix) 

in the sense of tempered distributions. 

Let X be a Brownian motion running with twice the speed of standard Brownian 
motion (i.e., X{t) ~ B(2t), where i? is a standard Brownian motion). The density 
of X{t) is -^e-l"=l'/4* (we denote by Kf^o,,^ = -^e'l^l'/^* the density of a 
one-dimensional BM starting at 0). From (|0.12p it follows that for /3 = ^ 

* 1 2i / 

= 2i5i/2(t/a;2)x-3 = exp 



x^^Ant^/x^ V 4t 

2 



(0.14) =^exp(-^ 



V 4t 



\B\ 
t:0, 



V2 

where KI..q\. is the transition density of \B\. Hence we get from the above observa- 
tions that 

(0.15) 2 r(2) 

'JtJ^t;0,x — '-'x^t;0,x- 

More generally, we have the following differential facts for K^.q^ (see equation (3.18) 
in [IS]) 
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Lemma 0.2 (Higher order PDEs associated with the A density). For (3 = 



E {2, 3, 4, . . the density K^^f^ satisfi' 



i^rr&^Kt^^^ t > 0, X > 

(-1)' 



(0.16) 



'-'x "^t-O.x 



(fc+1) 



x=0 

hm d'^K^f, 



.=0 r(i - 

= 0, t>0; 

0, t >0,k = 0,1,2,- ■■ 



t>0,k^0,l,2,--- ,{v-2) 



Let A^^\ . . . , A^") be n independent copies of a /3-inverse Levy motion of index 
< /3 < 1 starting at and independent of an n-parameter (n > 1) M'^-vahied Brow- 
nian sheet W^. We define the 7i-parameter R'^-valued inverse-Levy-time Brownian 
sheet (ISLTBS) based on and A^^', . . . , A^"), and starting at a; € W^, by 



(0.17) W^<„_. ._^(„,(t):=iy-(^AW(fi),...,A(")(t„)j; t € M!^. 

Clearly, P W^d, = ^ = 1 for t e dW^. Our first result in this more 

general setup gives an interacting time-fractional PDEs system that is solved by 
running the inverse-stable-Levy-time Brownian sheet in (|0.17p . It generalizes The- 
orem [U^] from the case /3 = ^ to that of < /3 < 1; and it also generalizes The- 
orem [Uini by giving other high (21/) order interacting PDEs systems for /? = l/f, 
u€ {2,3,4,...}. 



Theorem 0.3 (The ISLTBS fractional and high order systems). Let < /3 < 1. 

Let {W;^(i) S "■'^ n-parameter W^-valued LSLTBS based on 

a Brownian sheet and fi-inverse Levy motions {A^''}"_^ and starting on 
dWl at X e R'^ on {n,^,P). A ssume the same conditions of Theorem \0. 0\ on 
/ : R'^ ^ R. // 



u(t, a;) 



(0.18) 



r(^')(t,a;) = E; 




A(i)....,A(")> 



'A(i)....,A(") 



(t) 



j e N„, 



n ^^'Hu)\f 

^■»eN„\{i} / 

for {t,x) G M" X M'', then the family {u,'y^^^t, x)} ^^^^ is a solution to the 
fractional PDEs system 



(0.19) < 



'(a) a,% = iA,rW; 

(6) uit,x)^ fix); 

(c) rW(t,.T)=0; 

(d) rw(t,x) = [n,eN„\m^(/''i)^f] /(■^); 



, X e 



t e R 

t e dl 

tj = 0,.T e R'', 



i'l,x e 



for J e N„; where E{P, 7) = Ep(A(l)T) = j^^ x~^^gp{x)dx, for 7 > -1. 
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// we now assume 13 = l/u, v ^ {2, 3,4,.. .}, all the 21^ — 2 order derivatives 
of f are bounded and Holder continuous, with exponent < a < 1; (ind 



(0.20) '^}'\t,x) =Ep 



n A«(tO) /(wXa),...,Aw(t) 



V*6N„\{j} 



for (t,a-) e Wl X W^; then the family {u,'^y\t,x) \ is a solution to the 
iv-order PDEs system 



(0.21) < 



dt, ^ 2 



t g 



t/g dWl-\x<=, 



tj = 0, X- g M'', 



(6) u(t,x) = /(x); 
(c) '^J^')(t,.T) =0; 
[(d) '^J^')(t,x) = /(x)^^^'^(t); 
for j g N„; where 

ivw(t)= n Ep(A«(to)'=£^(i/i^,i'r"' n ^ 

ieN„\{i} ieN„\{i} 



(0.22) 



M(^)(t) 



»eN„\{j} 



.(A«(t,))^ 



n 

*eN„\{j} 



v/j ^ £^(1/^^,'^)' 



1=1 



and where ¥,p[A{t)'^] for — 1 < 7 is given by 

Ep[A(i)T] = ]Ep[(i(l)A)-T/''] = t'^/''Ep[L(l)-T/''] iT/''S(l/z/, 7) < 00. 



0.5. Conditional equivalence of the ISLTBS high order and time-fractional 
interacting PDEs systems. Our fourth main Theorem provides the conditions 
needed for equivalence to hold between the high order and the fractional systems 
in Theorem 10.31 



Notation 0.2. The fc-iterated /3-fractional derivative in t is denoted by 9, 



Theorem 0.4 (Conditional equivalence of the fractional and the (2i^)-order 
interacting PDEs systems). Assume that /3 = l/iy, v g {2,3,4,...}, all the 
2v—2 order derivatives of f are bounded and Holder continuous, with exponent 
< q; < 1 . Assume that, for every j g N,i , a function u-.W'ly.W^ ^M. satisfies 



(0.23) d\ 



r(^)£^(^.fc)"A^/(x)n.gM„\m^»^ 
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for fc = 1, • • • ,i> — 1. Suppose that, for every v, the pair {"^J , Y^^^ } satisfie 



(0.24) 



tj = 0, x e 

tj = 0, x € 



(a) 2''-^ A, 

(b) '^.^'\t,x) = f{x)Ni'\t); 

(c) rw(t,x) = /(x) n ^(^'i)^f; 

^{d) r(j')(t,2;) = ^J^'^(t,a;) = 0; t/eyK:;"%xe 
/or some j* G N„ (aH i/ie derivatives exist and are continuous) . Then (|0.19p is 
satisfied for j* by {u,'f^^">} iff ^filT) is satisfied for j* by ^u,'^y^Y ^'ori- 

versely, if fOM]) . fOTQl) . and ((0:2T1) are satisfied by {w, r^-^')} anrf |it,^j-'''| 
for j* e N„; t/ien (j0.24p holds for j* . In particular, under the condition (j0.23p . 
i/ie order system { (|0.2ip }^gpj and t/ie fractional system { (|0.19p }jgj.j are 
equivalent iff the collection of functions {"^j/i ^'■'^ j^gj^ satisfies (|0.24p /or 
j € N„. // {'^^^'^^^'Ijgf^ are t/ie ISLTBS functional defined as in Theo- 
remHa t/ien {^^J , ^(j)}^.^^^^ safe/y ([0:24| and ([0:23| /or ewer?/ j e N„. // 

n ~ 1, thenu — = 'V'^^'^ , the one-parameter version of the conditions (|0.24p 
and (j0.23p are satisfied and the equivalence between (|0.19p and (|0.2ip holds. 



0.6. Examples of the initial function /. Theorem 4.5 m Gilbarg and Trudinger 
[21] gives the fact that when /o is Holder continuous of order < a < 1 with 
compact support and / € C^(]R'') (with compact support) satisfy Poisson's equation 
A/ = /o, then / G C^'"(M'^) (the space of functions such that all the derivatives 
of order 2 are Holder continuous of order < a < 1 with compact support). 
Similarly, if /o € C^'"'{R'^) and / G C^{R'^) satisfy A/ = /o, then / e C^^+2'«(R'^). 
The condition of compact support can be removed by various means (see Theorem 
4.6 in Gilbarg and Trudinger [21] and the discussion before). 

We now give two related types of functions satisfying the conditions on the initial 
function / in our results. Let C be an arbitrary constant. Let /o M be given 

by 



(0.25) fo{x) 



then fo{x) is a-H61der continuous with compact support. When a = 1, /o G 
(joo -^yj^j]^ a,ll the derivatives being Lipschitz. We can then set / = /o in 
Theorem 10. 0[ Theorem 10. 1[ and Theorem 10.31 to get the corresponding PDEs. In 
the case of Theorem 10. 0[ if d = 2, then 



x\^ + \x 




4C(| 
A,/(x) = - 
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Alternatively, iterating the Gilbarg and Trudinger procedure above, we can use /o 
in (|0.25|) to generate functions whose derivatives of order 2v — 2 are all a-Holder 
continuous for ^ = ^ € {2, 3, . . .}. In particular, we can pick / such that A/ — /q 

and / e C^'"(R'^). In this case, with 



u(t,x) =Ep / {\ 



(t) 



(0.26) 



V»eN„\{j} / 



<.,(t) 



3 e N„, 



for (t,a;) £ M" x M'', we have that the family {xi.'^i/^^)^ .^^^ is a solution to the 
system of interacting fourth order linear PDEs: 



(0.27) 



(6) u(t,a:) = /(a;); 

(c) '^(J')(t,a;) 0; 

(d) '^«(t,z) = [n,eN„\m*'l /(^); 



t e 



t," = 0, X G 



for J G N„. We can also write the results in Theorem 10.31 in a similar way. 



1. Establishing the BTBS fourth-order and fractional PDEs systems 

LINKS 



1.1. A stochastic analytic unifying proof of Theorem 10 . l] and Theorem lO.Oi 

First, we need the following lemma that allows us to pull the Laplacian and bi- 
Laplacian in and out of different integrals. 

Notation 1.1. We bring the reader's attention to the difference between the spatial 
Laplacian 



(1.1) 

and the Laplacian 
(1.2) 

for sufficiently smooth i 



d 



A.g(u;(t,x)) 



d 

fe=l 



and g 



I, and for w : 



The second Laplacian (|1.2p appears when we apply Ito's formula in our proofs. 

Lemma 1.1 (Interchanging A, A^^, and A^ with different integrals: the BTBS 
case). Assume the same conditions on f as in Theorem \ 0.(A Then, recalling the 
notational conventions in Notation\0.1[ we have 



(a) Ep [a/ {w-^j. 



(b) A, 



Ep 



AlEp 



/ (l^(P.)) 
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(c) 

n 

i=l 

n 

and 

(d) 

n 
n 

1=1 

Remark 1.1. In parts (a) and (b) we wrote Wg^{pj) in place of since this is 
the form that appears in the main proof of Theorem 10.01 and Theorem 10.11 upon 
applying Ito's formula in (jl.7l) . making it more convenient for the reader to see the 
correspondence. 



^•4 


n •' . 




*eN„\{j} 



/ 


n 




ieN„\{j} 



/ 


n 




i6N„\{j} 



/ 


n 




*eN„\{j} 



Proof. Let if^^r"'^,!/ be the density of a Brownian sheet W^^{rj). Part (a) follows 
by integration by parts twice and the facts that 

hm dyJ{y)K^^i^.:^l^ hm f{y)dy,K^^j-:^l = Q 

which give 



(1.3) 



k=l 



k=l ' 

= A.Ep [/(M^4(r,))] , 
where the last equality follows from standard analysis and the fact that 



dy < oo 



for each k = 1, . . . , d by the boundedness of /. Part (b) is easily seen by noticing 
that, under the boundedness of the second partial derivatives of /, hsy{pj,x) = 

Ep A/ (w^^{pj)j is twice continuously differentiable in x; in fact, by the Brown- 
ian sheet-PDE connection in Lemma 1.1 in [T] hgy is the unique bounded classical 
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solution to the initial value heat equation with parametrized conduction coefficient 



(1.4) 



Thus, 



n 

»eN„\{j} 



MO,x) = A,/(x); 



dpj < oo, and hence 



A:, 



A/ ( W,-/p,) ) dp, = / A.Ep A/ ( W^^^ip,) ) dp 



' sy\H], 



which when combined with part (a) implies (b). Part (d) is exactly Lemma 2.1 in 
[Tj and part (c) follows the same steps with obvious modifications (see the more 
general case in Lemma [13] (a))- □ 



We now give a proof that simultaneously establishes Theorem 10.01 and Theo- 
rem [UTI capturing the fourth order PDE (|0.5p and the fractional PDE (|0.9|) in one 
equation. 

Proof of Theorem \0.0\ and Theorem \0.1l From our definition of the BTBS 
Worn Df„i we have that B'-^\ . . . ,5^") are n independent copies of a standard 

one-dimensional Brownian motion starting at 0, and we also have that B^^\ . . . , B*^"^ 
are independent of the n-parameter R''-valued Brownian sheet W^. Using these 
facts, we easily get that 



(a) u(t, x) = Ep / 

(b) '^(^■)(t,x) =Ei 



'b<i>,...,B(") 



:t) 



(1.5) 



(c) r'-j\t,x) = Ei 




Er[fiW^n]U^^Ld^' 

i=l 

f n /(w^<i),...3(">(t)) 

\*eN„\{j} / 

n 



n 

ieN„\{j} 



*eN„\{j} 



'B(i) B(") 



(t) 



n 

«eN„\{j} 



for j G N„. Now, fix an arbitrary j E N„, and let 



(1.6) 



Ew 



/(M^s>.)) 



0<p, <5,. 



We apply Ito's formula twice to / (W^) in (jl.Sp (a) in the j-th temporal parameter 
(see the detailed discussion in Appendix A in [T]); using the independence of the 

coordinates of the Brownian-sheet W^, the fact that |w^s^(pj);0 < < oo| is 
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a martingale with quadratic variation (w^^{-)'^ = njgN„\{j} ^'^"^ taking 

expectations and using the fact that the expectation of the stochastic integrals are 
zero, along with Lemma 1 1.1 1 (a) and (b), we get 



n 

*GN„\{j} 



(1.7) 



2" / 


n 




.»eN„\{j} 


2" / 


n 







S ■ 



n 

«eN,.\{j} 



A^fix) 



< 2" / 


n 


2 


> 




.*eN„\{j} 


Jo 


^0 «=1 



Equation p.7p contains in it the two BTBS-PDEs facts that {u, solves the 
fourth-order PDE in ([05]) and {m,^^-'')} solves the half-derivative PDE in 
We now show this fact. 

The interacting fourth order PDEs system in Theorem lO.Ol is seen by first taking 
the j-th time derivative dtj , with respect to tj , using the last equality in (|1.7p 
and the dominated convergence theorem to bring dt^ under the integral. We then 
use the fact that dtjK^^Q^. = ^d'^^Kf^^^. together with two integration by parts 
in which the boundary terms vanish. This last fact is because of the exponential 
decay of dg^ Kf^Q ^, and Kf^^Q ^, which nullify the polynomial growth of the inside 

integrals that is due to the continuity of A'^Vsy{pj , x) = A2;Ep A/ (^W^^{pj)j — 

which follows from Lemma 11.11 (a), the boundcdness of A/, and the discussion 
leading to the Brownian sheet PDE (|1.4p — and hence its boundedness on [0,rj]. 
We finish with two applications of the fundamental theorem of calculus, along with 
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Lemma [LT] (d). to get 

dt.u{t,x) = 



n 



f 


n 


J" I 






Jo Jo 



(1.8) = 



A./(x) 





n 'f 


n 




.jeN„\{j} 





./E™5a,/(x) + lA2'^«(t,x). 



On the other hand, to see the interacting fractional PDEs system of Theorem 10. 1[ 
we take the j-th time Caputo half-derivative d^. , with respect to tj , using the next 
to last equality in (|1.7p and the dominated convergence theorem to bring d^^. under 
the integral; we then use the fact that d^Kf.^^^^. = -^ds^KfJ^^ ,._ ^Soisjj^ 



(this was proved in equation (jO.lSp and the discussion preceding it) together with 
integration by parts, observing that J^^' AxVs^{pj,x)dpjSQ{sj)dsj = 0, along 
with Lemma ri. II (c). to get 



dlu[t,x) = -\T 





n 


1 JVJl 


jeN,.\{i} 



(1.9) 



= ^ ^ 2" 
\/8 



n ^» 

i6N„\{j} 



A^v,^{p,,x)dp,dl'WK^^l,A^ 

n 



^_LA,r(^)(t,x) 



The temporal-boundary conditions in (|0.5p and (|0.9p follow easily from the defini- 
tions of {u, , '^(■')} and the proof is complete since j is arbitrary. □ 



1.2. Proof of the conditional equivalence result: the BTBS case. Wc now 

establish the conditional equivalence result for the BTBS fourth order and fractional 
interacting systems. 

Proof of Theorem \ O.Si Fix an arbitrary j S N„ and suppose that two pairs 



'^So{sj) = linifjx^o ^t"^o s ^^^^ sense of tempered distributions {5Q{sj)dsj is the Dirac mea- 
sure at 0). 
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{uji^^^^] and {u,^^^] satisfy (|0.5p .and ()0.9p . Take the Fourier-Laplace transform 
{x k, and tj Sj) of (|0.9p to get 



(1.10) 



-\k\\,.u) 



Sy^Utyi-Sj, k) - S. ^^^f{k) - ^oj, 



Now, multiply both sides of equation (|1.10p by to get 



3 

1/2|, 12 



(1.11) 



,jUt^(Sj,fc) - /(fc) = 



V8 



Adding and subtracting .s~^/^ J ^'^t^^'l *' r{l/2){~\k\^f{k)) to (fTTTl) gives 



s,St^(s„ fc) - />) = ^^^(g,, k) 



(1.12) 



riieN 


Ab} 




24- 






riiGN 


Abl 





r(i/2)(-|fc|V"(fc)). 



Observing that 



lim A^y^^\t,x) 

ti—>-0 



(1.13) 



n i'l" 



^eN„\{i} 



A,/(.T 



2" ^n»eN„\{j}^' 



r(l/2)A,/(a:);a:e 



equation (|1.12p inverts to 
du 



(1.14) 



at,- 



y^P^A./W + ^A.(5i/V(^)(t,.));t e e R'^. 



Comparing equation (|1.14p and equation (jO.Sp (a) we must have 

-^A,(a^/V(^'(t,x)) = lA2'^(^)(t,x);t eRl,xGR'' 

and condition (|0.1ip is satisfied. 

For the rest of the proof, suppose that (jO.lip holds. We then assume first that 
{u, yii)] is a solution to p^Q]) . Observe that 



dli'^{dli'^u{t,x)) = dt,u{t,x) 



%^A./(x) 



This follows from the fact that, under the smoothness conditions on u, it satisfies 



-1/2 



dli'^i.dl'^'^u{t,x)) = dt^u{t,x) 
and the fact that 

<\(t,.)|,.o = ^lhn„A.yO-)(t,.) = -l 



dl'^'^u{t,x)\t^^o 



1 / n.gN„\m^' A X 
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Taking the half derivative with respect to tj of the equation in (|0.9p then we get 
that {u, ^■'} is a solution to (jO.Sp . 

If we now assume that {u,'^^] is a solution to (|0.5p , then we can use condition 
(|OlT|) (a) to rewrite (|03|) as 



Taking the Fourier-Laplace transform (a; k, and tj Sj) of p.lSp using (|0.1ip 
(d) we get 

(1.16) 



s,ut,{sj,k)~f{k) - i/npiJM^.7V2r(i/2)(_|fc|2/(fc)) 



Observe that the first and third terms on the right hand side of (|1.16p cancel each 
other. Hence (|1.16p reduces to 

(1.17) s,u,^(^s,,k)- m = ^Hk\')sf%f{s,,k). 

1/2 

Dividing both sides by Sj we obtain 

sfut,is,,k)~ s-'/'m ^ ^Hkf)%f{s,,k)^ 
which, upon inversion, implies that {u, '^t-')} is a solution to (10.91) . □ 



2. Establishing the ISLTBS high-order and fractional PDEs systems 

LINKS 

2.1. A stochastic analytic proof of Theorem 10.31 Like the BTBS case, we 
need a lemma that allows us to pull the Laplacian and bi-Laplacian in and out of 
different integrals. Let A(i) be the inverse of L{t) with index < (3 < 1. We can 

calculate Ep[A(i)'^] for -1< 7 as follows: A{t) = (L(l)/t)-i/'^ by Corollary 3.1 in 
Meerschacrt and Scheffler [55] hence 

Er[A{ty] - Ep[(L(l)A)-^^] = r^Ep[L(l)-^^] ^^£{(3,^) < 00. 

Lemma 2.1 (Interchanging A^: and with different integrals: the ISLTBS case). 

Let f : R'^ — > M &e bounded and measurable such that 
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(a) / bounded and Holder continuous, with exponent < a < 1, then we have 



A, 



n 

*eN„\{j} 



n 

»eN„\{j} 



n 



1=1 



(b) all the 1v — 1 order derivatives of f are bounded and Holder continuous, 
with exponent < a < 1, then we have 



A" 



n 

jeN„\{j} 



n 

*eNn\{j} 



n 



Proof. The proof of (a) is similar to the proof of Lemma 2.1 in [I]. We need to 
consider second derivative instead of fourth derivative and the function / instead 
of the second derivative of h there. We have 



(2.1) 



dlEr[fiW^n]= f{y)dlKf.^l}--'')dy 



(xfe-iy-(s))2-n:=i5. 



n:=i. 



{nwn~f{x)) 



where Kf-^x^y^''^'' is the BS density previously given in (|0.2|) . By scaling we obtain 



2 n 



sjEf 



1=1 



for some constant C . Using (|2.ip and Cauchy-Schwarz inequality, along with the 
fact that / is Holder continuous with exponent a, we get 



(2.2) 



diE^ [f{w:)] 



< 



c 



nil 



l-a/2 
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Now, 



< 



n 

*eN„\{j} 



n 

ieN„\{j} 



i=l 



c 



= 1 •^i i=l 



(2.3) 



C 



n ^ 

*eN„\{j} 



a/2 



,-{l-a/2) 



/ /"CXD \ poo 

= CEp[A(t,)"/''^] n Ep[A(tO^"/')] 
*eN„\{j} 

= Ct5"/'"'^/"^(l/i^,a/2 - 1) [| tf/('")£;(l/zy,a/2) < oo 

jeN„\{j} 

We next prove part (b). For notational simplicity we show that 



(2.4) 



n -^^ 

'6N„\{j} 



n 

*6N„\{j} 



n 



the mixed partial derivatives cases can be proved similarly. 
Using the boundedness of / and by observing the facts 

^^hm^ di;'f{y)c}^:-^K^^;^y = 0, J = 1, 2, . . . , (2Z. - 2) 

(since all the derivatives up to order 2i/ — 3 are bounded and Lipschitz in yk) and 
using integration by parts (2i/ — 2)-times we get (using symmetry of K^-^x^y''^'^ in x 
and y) 



(2.5) 



9^rEr [f{W^)] = I f{y)dl:Kf.,l^--^^dy 
f{y)dl:Kf.f,}-'''Uy 
dlr'f{y)dlKf^!^^'^dy 

(xfe-M/-(s))2-n:u«. 



{h{W^)-h{x)) 
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Where 9,^^ ^.fiu) = h{y) is Holder continuous with exponent < a < 1. Hence we 
can finish the proof of part (b) as the last lines of the proof of part (a). Now 



< 



n 



n 

«eN„\{j} 



dtl¥.^[f{W:)] 



i=l 



c 



= c 

(2.6) 



= C 



n - 

*eN„\{j} 



u-l+a/2 



I-q/2 



i=l 



^(l-a/2) 



1=1 



n 



i^-1+q/2 7^A,/3 



= CEp[A(tj)-(i-"/2)] ]J Ep[A(t,)*''"^+"/^^] 
'eN„\{j} 



«eN„\{j} 

hence we obtain (|2.4p by a standard classical argument. □ 



< oo. 



Theorem 10.31 can now be proved. 



Proof of Theorem \0.S[ Throughout the proof j € N„ is fixed but arbitrary. Let 
< /3 < 1. Let Vsy{pj,x) be given by (|1.6p . Following the same steps of our proof 
of Theorem 10.11 above gives 



u{t,x)=Ep [/(w^(i),...,A(")(t) 



(2.7) 





n ^» 




.ieN„\{i} 



Er[fiw,n]YlKisA^ 



wc then take the j-th time Caputo /3-derivative 9^. , using the dominated conver- 
gence theorem to bring under the integral; we then use the fact that 9^. K^,'.q = 

—dsjK^.'^^^ — yIi-p) '^o(gj jl together with integration by parts, observing that 



■^(5o(sj) = linifjx^o ^t'<) s ™ ^^'^ sense of tempered distributions (5o{sj)dsj is the D 
sure at 0). 



irac mea- 
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■fo° -fo^ ^xVsj,{pj,x)dpjSQ{sj)dsj = 0, along with the fundamental theorem of cal- 
culus and Lemma Em (a), to get 



{a 


n ^» 




.*eN„\{j} 



(2.8) 




n 

*6N„\{j} 



1=1 

n 



i=l 



The temporal-boundary conditions in (j0.19p follow easily from the definitions of 
{u, Y^^")}, and the proof of our ISLTBS connection to the fractional PDEs system 
in (|0.19p is complete since j is arbitrary. 

Now, fix an arbitrary ^ \/v such that v E {2, 3, 4, . . .}. Define the simplex 

§j {(n,T2, ...,T,) e [0,s,]'';0 < ri < r2 < ... < < s^} . 

We use the notation 



k=l 



h{Ti)dTi ■ ■ ■ dTiy 



Let Vsy{Ti,x) := Ep f [W^^{Ti)j (the j-th temporal parameter is ti). Applying 

Ito's formula ly times to f (W^) in (|0.4p . in the j-th temporal parameter, using 
the independence of the coordinates of the Brownian-sheet , the fact that the 

quadratic variation (W^^{-)) = t'j HieN \{j} ^^'^ taking expectations, along 

with repeated use of Lemma [TTT] (a) and (b), we get 



it,x)^f{x) + J2 



A^/(x) 



(2.9) 



fe=i 







h 


n 




.jeN„\{j} 



k n 





n < 




.'eN„\{j} 



Kv.M^x)<^dTkX{Kt:.l,As 



fe=i 1=1 



Taking the j-th time derivative dt^ , with respect to tj , in (|2.9p and the domi- 
nated convergence theorem to bring dt under the integral over Rlf. ; and then using 
Lemma 10.21 together with i/ integration by parts (in which the boundary terms 
vanish as in the BTBS case — using similar argument as the one before (|1.8|) in the 



BTBS case above, A'^Vsy{Ti,x) = A^E^ A"'^ f {w^^{ti)^ is continuous by the 

conditions on / and its derivatives of orders up to 2i/ — 2 and the conclusion follows) 
and V applications of the fundamental theorem of calculus, along with Lemma [2^1 
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(b), to get 



(2.10) 



k=i I ■^"+ \ieN„\{j} / ■ 1=1 



The temporal-boundary conditions in (|0.2ip follow easily from the definitions of 
{u, '^J}, and we are done since j is arbitrary. The terms Ep[A(t)'''] for — 1 < 7 as 

follows: First A{t) = {L{l)/t)-^ by Corollary 3.1 Meerschaert and SchefSer (2004) 
hence 



Hence 



.[Am = Ep[(L(l)A)-^'5] = t^''Ep[L(l)-^^] =: i^'';7(/?,7) < 

Ep(A«(t,))" 



MW(t) 



n Ep(AW(t, 

ieN„\{j} 



n 

*eN„\{j} 



1=1 



(2.11) 



uii/i^,Ky 



n 



»eN„\{j} 



n tC/(l/z^,i')= (Ep[L(1)-1] 

»6N„\{j} 



n 

»6N„\{j} 



completing the proof. □ 



2.2. Proof of the conditional equivalence result: the ISLTBS case. The 

following will be useful in the proof of Theorem 10.41 

Lemma 2.2. Let < /32 < 1 and /?i + /?2 < 1 then we have 



fil+132 



dtf{t) 



r(i-/3i) 

if all the derivatives exist and are continuous. If (3 — l/v, v G {2, 3, 4, . . .}, then 



(2.12) = dtfit) fTTT^ 

K=l ^ v' 

if all the derivatives exist and are continuous. 



fit) 



t=0- 
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Proof. This follows by a simple application of the Laplace transform of the Caputo 
fractional derivative. □ 



Remark 2.1. We can use Lemma [2.21 and the method in the appendix to show 
that the system {u,'^J^ in Theorem 10.31 satisfy (|0.21|) . First note that 



(2.13) 



2^- 



9fc 



n Ep(AW(t,) 

«eN„\{j} 



E{i/iy,kr~' n t: 

i6N„\{i} 



k/i^ 



Hence we can show that 



(2.14) 



^ ' .tt r(i - - k)/v) 



Of. " u 









^A^%^^^ + y -j J- 



AV(x) 
2fc 



ii;(i/^,fc)"-i n < 

*eN„\{i} 



Now comparing the equation (|0.2ip and (|2.14p we see that 



£;(l/i.,fc) = Ep(A(l)^-) 



iy{k-l)\ 



Proof of Theorem \0.4\ Fix an arbitrary j € N„ and suppose that two pairs 
and satisfy (|03T|) and (|0J9)) for /3 = Applying ^(pi)®!/'' 

to both sides of (|0.19p and using condition (|0.23p and Lemma give 



u-l ^-(v-k)/v 

^ ^ F(l - (z. - /c)/:/) 



o^® 1/ 









(2.15) 



i^2'=(fc- 1)! 



^ l^.^^O-) ^^^^-(.-fe)/.i?(l/^,fc)"Al/(x) 
^ fe=i 

- V 

^£^F(l-(i.-A:)H 



n 

'eN„\{j} 



T^A^^J-'') 



establishing (|a24)) (a) and thus (|0^ since (|0^ (b), (c), and (d) are satisfied by 
assumption. 

For the rest of the proof, suppose that (|0.23p and (|0.24p hold. We then assume 
first that {u, r^^)} is a solution to for (3^1/ v. Applying £)(j-i)®V>' ^^^^i^ 
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sides of (|0.19|) as in p.lSp and using the conditions in ()0.23p and (|0.24|) we get the 
fact that {u, '^^} is a solution to (laMj) . 

Now, let A^{t,x) := ^A^^'^'-') (t, x). If we assume that \^u,'^J} is a solution to 
dnSIl), then we can use (((1^ and ((DTM)) to rewrite (^I^ as 



(2.16) 



fe=i 



*6N„\{j} 



,-(i/-fe)/i/ 



nfc® 1/ z-* 

''J 









Using Lemma [2.21 we get 

(2.17) a.f'/'^M = 

Now, taking the Fourier-Laplace transform (a; — ^ fc, and tj — > s^) yields 



(2.18) sy''u,^{s„k) ~ sy^-'fik) ^ Ai^{s,,k) 

Taking inverse Fourier-Laplace transform implies that {w, '^'^■'^} is a solution to 
(|0T9l) for /3 = l/v. □ 



3. On RELAXING THE BOUNDEDNESS CONDITION ON / AND ITS DERIVATIVES 



As with some problems of applied mathematics, the boundedness assumption 
on / and its first few derivatives can sometimes be too restrictive. We now briefly 
discuss how to relax this boundedness condition. Carefully examining the proofs of 
the main results reveal what types of unbounded / are allowed. We need Lemma ll.ll 
to hold for the BTBS PDEs results to hold, and we need Lemma 11.11 (a) and (b) 
and Lemma O for the ISLTBS PDEs resuhs to hold. 

In addition to the Holder continuity on / and all of its 2;^ — 2 order deriva- 
tives, V = ^ {2,3, . . .}, we need some growth and integrability conditions with 



respect to the BS density given by (|0.2p . Namely, we can replace the boundedness 
assumptions on / and its first 2;^ — 2 derivatives by 



(3.1) 



(i) 



(ii) 



dy < oo; i = 0, 1, . . . , 2i^ - 2, j = 0, 1, . 

and i + j < with k,l — I, 



,2i/, 



Y.dlfiy)K?,'}:-:^ldy 



k=l 



dpj < oo; Z = 1, . 



, , d, rj > 0. 



Of course, in particular, p.ip implies 
(3.2) hm di-'f{y)c^:-^K^ 



n,d) 



0; ] - 1, 



,(2z.-2), fc = 1, 



By the exponential decay in the BS density and its spatial derivatives, condition 
(|3.1[) holds when / and all of its derivatives of order up to 2;/ — 2 have polynomial 
growth. 

The interested reader can check that Lemma [TTT] and Lemma [2TT] hold if we pre- 
serve all the differentiability and Holder continuity conditions on /, while replacing 
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the boundedness on / and its derivatives by p.ip . This means that, without sig- 
nificant changes in our proofs above (except obviously using condition p.ip on / 
and its derivatives in place of the boundedness condition on / and its derivatives 
of order up to 2u — 2), we have stronger versions of all our main results that we 
capture in the following theorem. 



Theorem 3.1. Assume that (3 = l/v, v e {2,3,4,...}. Then, Theorem 
Theorem \0.4\ all hold with the boundedness conditions on f and its derivatives 
of order up to 2v — 2 replaced by the conditions in (|3.ip . leaving all other 
conditions in force. 



Wc note that the BTBS results (Theorem [021 Theorem UTTl and Theorem 
fall under the case oi v = 1/2, with a minor scaling of the density {At is replaced 
by 2t), as in equation (|0.14p . 

Acknowledgement. We are sincerely grateful for the anonymous referee's careful 
and constructive comments which improved the readability of the paper. In partic- 
ular. Subsection 10.61 and Section [3] were added as an answer to an excellent point 
made by the referee. 



Appendix A. Another proof of the fractional PDEs connection in 

Theorem [QTsI 



Here, we present a different proof of the fractional PDEs connection Theorem l0.3l 
that is an adaptation of the proof of Theorem 1.1 in [1]. Take the time fractional 
derivative and put it inside the integral by the dominated convergence theorem, 
using part (e) of Lemma 10.11 and integration by parts once, and the boundary 

conditions lims^._j.oo K^.\q ^^-^ lims^^o -^'t^'o — ~ r(i-,3) ' ^'^'^ using Lemma 

1.1 in [1] to get 

diu[t,x)=j EAfiW^md^.^K^;^,^^^ n ^L'^ 

^^^Ep[/(M/-(s))] Ld.,,K^;l,^^^^^^5,{s,)\ n K^ls.d^ 



ieN„\{i} 



i=l 



(Ai) I r' \ 



*eN,.\{j} 



r(l-/3)^ 
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To justify taking the Laplacian outside of the integral we use Lemma |2. II Next we 
verify the boundary condition (d), since at tj = 0, Ep [/ (Ty^(s))] = f{x) we get 



(A.2) 



= [ A li A /(^) n ^tLd^ 

■^^+ \ieN„\{j} J ieN„\{j} 

= fix) n (A(to) 

«eN„\{j} 

=/(^) n (c(/^)^0 

«eN„\{j} 

Where the last equality follows from Corollary 3.1 in [55] 



Appendix B. Frequent acronyms and notations key 
I. Acronyms 

(i) BTBM: Brownian-time Brownian motion. 

(ii) BTBS; Brownian-time Brownian sheet. 

(iii) BTP: Brownian-time process. 

(iv) ISLTBS: inverse-stable-Levy-time Brownian sheet 

IL Notations 

(i) N„ = {!,..., n}. 

(ii) Kf.^Q 3 . : The density of a one-dimensional BM starting at (see Section 
Ell'just before (pij)) ). 

(iii) K^.'.Q J. . : The density of a one-dimensional /3- inverse Levy motion A 
astarting at (see Section ing . 

(iv) K^.^^y''''^ : The density (or kernel) of an n-parameter d-dimcnsional 
Brownian sheet (see (|0.2p ). 

(v) Ax vs. A (see Remark I l.ip 
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